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ABSTRACT 

Xncooq^essible txirbiiLent flows are investigated in the frame- 
work of ideal amgnetohydrodynamics. All the field quantities vary 
with only two spatial dimensions* Equilibrium canonical distributions 
are determined in a phase space whose coordinates are the real and 
imaginary parts of the Fourier coefficients for the field variables. 

In the geometry considered, the magnetic field and fluid velocity have 
variable x and y cooponents, and all field quantities are independent 
of z. Three ccmstants of the motion are found (one of them new) which 
survive the truncation in Fourier space and permit the construction 
of canonical distributions with three independent ten^eratures. 
Spectral densities are csdculated. One of the more novel physical 
effects is the aj^earance of macroscopic structures Involving long- 
wavelength, self -generated, magnetic fields ( '^magnetic islands") for 
a wide runge of initial parameters. Current filaments show a tendency 
toward consolidation in much the same way that vorticlty filaments 
do in the guiding-center plasma case. In the presence of finite 
dissipation, energy cascades to hl^er wave nximbers can be acccmipanied 
by vector potential cascades to lower wave numbers, in much the same 
way that in the fluid dynamic (Navier-Stokes) case, energy cascades 
to lower wave numbers acojmpany enstrophy cascades to higher wave 
numbers. It is siiggested that the techniques may be relevant to 
theories of the magnetic dynamo problem and to the generation of 
megagauss magnetic fields wl:^n pellets are irradiated by lasers. 
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1. INTRODUCTION 

We have recently been Involved in a miniher of investigations 
in the statistical theory of txirbulence for two-dimensional flows 
in electrostatic guiding-center plasmas and/or inviscid Navier-Stokes 
fluids. Any tractable mathematical description of a turbulent con- 
tinuum requires a discretization^ and this has been carried out in 
one of two ways; either in terms of long rod-like "particles" 

(Onsager 19^9; Taylor and McNamara 1971; Montgomery 1972; Jc^rce and 
Montgomery 1975; Montgomery and Joyce 197^; Seyler 1974; Montgomery 
1975a, b; Lundgren and Pointin 1975; Pointin and Lundgren 1975), or 
in terms of truncated Fourier series representations of continuum 
fluid equations (Kraichnan 19^7, 1975; Seyler, Salu, Montgomery, and 
Khorr 1975; Montgon»ry 1975a, b; Montgomery and Sedu 1975)* 

The two discretizations lead to dynamical systems which differ 
in their predictions at the anallest spatial scales, but the most 
interesting physical phencaaena have occurred at large spatial scales 
and have been comnran to both representations : for a wide range of 

initial conditions, the flows have shown a remarkable capacity to 
organize themselves into a few large, persistent vortices, comparable 
in cross-sectional area to the dimensions of the system. In the 
plasma interpretation, phase space considerations dictate that for 
the higher-energy states of the system, the more probable configurations 


4 


consist of large regicms of non-zero net charge density, around which 
the plasma E X B drifts azlmuthally in a direction nomal to the dc mag- 
netic field* These configurations have repeatedly appeared in ntimeri- 
cal simulations (see the above references) and have also been studied 
from the point of view of time -independent self-consistent Vlasov 
equilibria (Book, McDonald, and Fisher 1975 J McDonald 1974). 

The techniques by means of which these phenomena have been 
investigated are of considerable generality, and the limits of their 
applicability have not been reached. Here we report sc«ne theoretical 
investigations into a situation in which the plasma flow is still 
two-dimensional, but with a hi^ enough ratio of mechanical energy 
to magnetic energy that the magnetic field has to be self-consistently 
determined: colloquially, the "low beta" limit is no longer appro- 

priate. (Heretofore, the plasma has been assumed to E x B drift 
across a given, static magnetic field, normal to the plane of varia- 
tion of the field quantities.) Some important beginnings on the 
problem have been made by Frisch, Pouquet, Leorat, and Mazure (1975) 
who considered the three -dimen sicaial magnetohydrodynamic case, with 
conclusions significantly different from the ones we report here for 
two dimensions. A two-dimensional calculation with spatial variation 
across an external dc magnetic field is due to Schumann (1975). 

We proceed by means of the equations of ideal magnet ohydrod 5 ma- 
mics, and consider incompressible flows in the following geometry. 

There is no externally Imposed component of the magnetic field, but 
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there may he a lar^ internally-generated one. !I%e ma^ietic field B 
and the fluid ve3x>city v have only x and y con^caients. All quantities 
are Independent of the z-coordinate. The electric current vector 4 
and the vorticity vector where uj = 7 X v and 4 - c7 X B/^, have 
only z-conponents. S(»Qe numerical simulations for this geometry have 
be«i carried out by Tappert (l97l)» At any stage in the development, we 
may add a uniform dc magnetic field in the z-direction without changing 
any of the manipulations. We omit it from the development to keep the 
equations as siirple as possible. 

The omission of dissipative effects (finite conductivity or 
viscosity, for exan^le) will be expected to alter the quantitative 
conclusions significantly, just as it does in the Navier-Stokes theory. 
It is c<xnmon knowledge that omission of viscosity, however small, 
alters the inertial range spectral behavior with k for the Navier- 
Stokes equation. [Compare, for example, the results of Seyler elb al. 
(1975) or Basdevant and Sadourny (1975) with those of Herring, Orszag, 
Kraichnan, and Fox (1974)*] But qualitative conclusions, independent 
of the presence of small dissipation at high wave numbers, are easier 
to draw in the non-dissipative limit, based as they are on the exis- 
tence of constants of the motion which decay slowly in the presence of 
a finite dissipation. The purpose of this article is not quantitative 
ccxnparison with data, which would in any case be premature, but the 
Isolation of some qualitatively new gross physical effects which have 
not to our knowledge been previously calculated. 

In section 2, the Lundquist equations for the present geometry 
are written down in an appropriate set of dimensionless variables. 
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Each field quantity is represented by a Fourier series involving a 
large but finite number of terms. The Fourier <»efficients are then 
the basic time -dependent dynamical variables of the theory. Qaestions 
of the convergence of the Fourier series results to solutions of the 
original Lundquist equations are non-trivial, but they are as intract- 
able here as they have been in fluid turbulence theory and are left 
open^ as they must be at present. 

We work without reference to nonaal modes, eigenfrequencles, 
waves, perturbation theory, or any of the paraphernalia which have so 
unfortunately limited many plasma turbulence calculations to minor 
corrections to linear theory. There is no attempt to follow the 
evolution of the dynamics from an initial to a final state. Rather, 
the procedures followed are basically those of Kraichnan (1967) 

[see also: Seyler et al. (1975)]# in which classical equilibrium 

ensembles are built aroiu^ constants of the motion identified from 
the Fourier transformed equations of motion. Constants of the motion 
which exist for the original equations but not for their truncated 
Fourier representations are regarded as insufficiently rugged to 
enter into the statistical formulation of the problem. Identification 
of constants of the motion is at first glance more an art than a 
science, and whether one has found them all is always slightly uncer- 
tain. Nevertheless, symmetry ccaislderations, ccanbined with an accumu- 
lation of practical experience, lead to confidence that the possi- 
bilities have been exhausted after awhile. The striking confirmation. 
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by Seyler et etl. (1975)# of the Kralcbnan theory's predictions lend 
considerable credibility to the conject'ure that cmly tvro such constants 
(energy and enstrophy) exist for the truncated Ifavier-Stokes system. 
Until mathematical investigations of a ccaislderably higher order of 
rigor than those that have been carried out so far are performed, 
this is probably the best we can do. Similar reservations apply to 
the work of Frisch e^ al. (1975)* 

Once the constants of the moticxi represented in terms of the 
Fourier coefficients are identified, the statistical formulation of 
the problem is presented in a phase space whose coordinates are the 
real and imaginary parts of the Fourier coefficients. Canonical en- 
sembles are constructed in this phase space by classical argument. s 
(e.g., ter Haar 1967) in which the probability of a state is maxi- 
mized subject to the constraints that the constants of the motion 
retain their initial values. The method of Lagrange multipliers Is 
utilized, with the Lagrange multipliers playing the role of inverse 
temperatures, one for each independent constant of the motion. The 
derivation of the ensemble by this conventional Boltzmann-Gibbs 
combinatorial method differs in no important particular from the 
primitive derivation of the Gibbs distribution offered in any first 
course in statistical mechanics. From the canonical ensemble so 
derived, expectation values of any functions of the amplitudes of the 
turbulent field can be computed by straightforward integration (e.g., 
the spectral densities). In section 5 some of these are computed, 
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and their variation with initial values of the constants of the 
nation is discussed. 

Of particular interest to us has been the question of the 
degree to which a highly ojnducting turbulent fluid in which there 
is an initially small amount of magnetic field energy can convert a 
significant fraction of its turbulent mechanical energy into magnetic 
field energy— ‘l.e. , the degree to which the plasma can "magnetize" 
itself. In a number of situations (e.g., in the appearance of mega- 
gausn magnetic fields in the irradiation of pellets by lasers, or in 
dynamo the<n:les of the earth's magnetic field), spontaneous large 
macroscopic magnetic fields are known to appear in magnetohydrodynamic 
fluids. While our theory has not advanced to the point where a serious 
qviantltative confrontation with data can be undertaken for either 
of these phenomena, we do show in principle how such large internally- 
generated B fields can develop, a result believed to be among the 
first of its kind insofar an it is rigorously derived, albeit within 
a highly sinq)lified framework. 

A second phenomenon which can be extracted is the consolidation 
of electric current filaments, much in parallel to the consolidation 
of vortices in the guiding -center plasma in two dimensions. 

These and other results are ccsnmented upon further in section 
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2. BASIC EQUATIONS AND CONSTANTS OF THE MOTIOI 


The lundquist equations of ideal inconqpressible magnetohydx^- 
dynamics are: 


^ = 7 X (v X B) 
dt 


* 


\ 4 

■ c * 5 • ' 

7 • V = 0 . 

where = (O, 0, J^) is the electric current density, B = (b , B , O) 

X y 

is the magnetic field, v = (v , v , O) is the fluid velocity, p Is 

~ X y 

the pressure, p is the (constant) mass density, and c is the speed 
of lig^t. The vorticity uj = (O, 0, m ) is given by 




= 7 X V 


We assiune throughout that = 0, and that periodic boxmdary condi- 
tions apply in x and y. In the given ge<^etry, a uniform, constant. 
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magnetic field in the z-direction can be added without altering the 
equations of motioi. It simplifies the expressions to omit it| so 
we do. It can be re-inserted at any step of the de^lops^nt. 

We measure velocities in units of a constant velocity 
which may be taten to characterize the mean Initial speed of the 
turbulent field. We measure lengths in units of L^, which may be 
taken as a characteristic maczoscopic length of the turbulent field. 

measure tln»s In units of measure ma^etlc fields in 

units of s where is the mass density. We measure 

current densities in units of s cB^4ttL^» In these units, Eqs. 
(l) reduce to 


95 


(2a) 


9v 2 

— + v»Vv = (7xB)xB-cf7p , 


("’b) 


7 • V = 0 , 


(nc) 


j = 7 X B 


0) = 7 X V 


(T’d) 


2 ^ 

Here, c“ = p/aJJt is the ratio of the characteristic mechanical 
s o '^o o 

o 

pressvire p^ to As usual, taking the divergence of Eq. (?b) 
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and using Sq. (2c) enables us to solve a Poisscsi equation for p, 
so that p is in effect a knoim functional of B and v. 

All quantities are exp»»3ed as Fourier series in a lu'ge 
square box, assximing periodic boundary (Editions: 


B = E B(k, t) exp (ik • x) , 
k 


X = E t) exp (ik • x) , 

Is 

‘"z “ E • 3£) 

k 

^ k 


( 5 ) 


where k has only x and y cooponents; 
ik X e 

U) ■ J(i) , 

\ 

ik X e 

v(k) = «^k) , (U) 


and e_ is a unit vector in the z-direction. 
z 
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The dynamical equations can he written entirely in terms of 
the scalars t) and J(k, t), the Fourier transformed vorticity 
current density and electrical current density. A straightforward 
substitution of Eqs. (5) ard (U) into (2) yields 

dd^k; t) 

“it — “ (5) 


and 


dj(k, t) 

(C) 


6(j + r-k) = 1 when its argument is zero and is zero otherwise, end 


?) = V®' £) " I p) 


2 “ 2 

P r 


(7) 


P) " -Mp* r) = I Sz * X P) 


I 2 \ 
‘ k ' 


O o 

\r‘ p" 


( 8 ) 


The sums in Eqs. (5) and (6) are over all non-zero wave vectors p and 
r. Dropping Eq. (6) and the j-terms in Eq. (5) reduces the system 
to the inviscid Navier-Stokes equation in the same geometry (iJeyler et 
al. 1975). C^ly wave numbers peraitted by the periodic boundary ccnJi- 
ticxis are, of course, allowed. 
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Equations (?) and (6) desoribe an incooq^ressible flow in the 
Idiase space whose coordinates are defined by the real and imaginary 
parts of the u)(k, t) and the J(k, t). Hence a Liouville equation 
obtains in this j^ase space; we may seek equilibrium enseiA>le3 in 
this phase space in the conventional fashion, providing that we are 
able to find the isolating integrals (which define hypersurfaces in 
the phase spa<%) or "c<mstants of the motion". In particuj.\r. we 
are interested in the constants which stirvive the truncation which 
represents Eqs. (5), (?), and (6) in terms of a large but finite 
nu-c^r of terms* 


Many constants of the motion exist for Eqs. (.'?) which do not 
survive the restriction, in Eqs. (?) and ( 6 ), of values of |k(, |p| , 
and |r| to a large but finite maximum. In particular, the line 
integral J ^ B x dt, along any curve whose -nd points are 1 and :■ 
and which moves with the fluid velocity, is time-independent according 
to Eqs. (2) but not according to Eqs. (?) and (6), if any finite k 

rnsix 

is allowed to restrict the values of |k| , |p|, and jr]. 

We have found three constants of the motion for Eqs. ( r ) which 

do survive the restriction that |k| , |pj , |r| are required to lie in 

an annulus betv/een k . and k (k . is just 2 tt divided by the l>ox 

min max min 

size). We believe that they are the only ones, but cannot rigorously 
prove it. They arc the energy , 


e 


4 /(/ + B") dxdy 


{ ‘’5a) 
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tl» "dross helicity" P (Woltjer 1958, and Prisch et 1975), 


P H 



B • V dxdy 

rv/ /Srf 


(9b) 


and a new integral, the mean sqtiare vector potential. 


^ = > (9c) 

where B = ^ X and ^ is the vector potential in the gauge for 
which 7 • ^ = 0. The integrals in (9) rxin over the basic box. 

The Fourier series representations of g., P, and A are 


e = 





r 


P = 


1 

2 


E 


4k) 3(-k) 


> 



( 10 ) 


(The time arguments are conveniently emitted.) The sums in Eqs. (lO) 
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are over all the k values allowed by the periodic boundary conditions 
between k . and a large but finite k„^^, which may be chosen arbi- 
trarily. It is an easy manipulation to show that the expressions 
given in Eqs. (lO) are conserved by Eqs. ( 5 ) and ( 6 ). 

The canonical distribution, to be inferred by entirely con- 
ventional arguments, is then 

D = T] exp (-ae - pp - rA) , (11) 

where D is the probability distribution in the many-dimensional phase 
space whose coordinates are the real and imaginary parts of 

‘"(is) = ‘"i.(k) + 

and 

J(is) = jp^is) i^i(ls) • 

T] is a normalizing constant to be determined by the normalization 

where f dX is an integral over all the independent oi , tv, , 1 , 1 . 

r 1 r * i 
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[Since <u(k) = u) (-k), j(k) = j (-k). Independent are 

associated with (»ily half the k-vectors. ] a, S, 7 play the roles 
of inverse ten 5 >eratures, and are constrained by the requironents 
that Eq. (ll) be normalizeable— i.e., that for large | j(k)( ' and 
lu)(k)t^, the argument of the exponential is a nK>notonically decreasing 
quantity. This implies right away that 

a > 0 . (15) 


Both 7 and ^ may be negative, but not a. 

Another class of constants of the motion which do not survive 
the truncation of the Fourier series are integrals of the form 

= J (7^ dxdy, where n is some integer and the integral is over the 
basic box. The constancy of can be proved by using the governing 
equation for the vector potential. 


Jl 

at 


+ y • 



However, only for n = 2 does remain constant in the face of the 


truncation. 
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5. CHARACTERISTICS OF THE DIFFERING 
!miP£RATUR£ REGIMES 

S,, Pf and A are stuns of terms, each of which is indexed by a 
single wave number; therefore Eq. (ll) predicts that in equilibrium 
there is no correlation between coefficients corresponding to dif- 
ferent k’s. D factors into a product of distributions, one for each 
k, and the single-k distribution function is 

Jj.(k), Jj^(k)] 

/V 

= exp + cu^(k) /fc2 

+ o^(k)]/k^j , (lU) 

where Tlj^ is a normalizing constant, given by: 
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Cie additional Inequality required in oi^er that the distri- 
butions be normalizable is that 


a > 




(15) 


for edJL k. If 7 < 0, Eq. (l5) requires that 


a > 




{ l6a) 


while if 7 > 0, Eq. (15) requires that 


a > 




(l6b) 


Ejcpectation values computed with respect to Eq. (l4) are: 


<J®(k)> = (3f(k)> = 





(17) 


and 
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(18) 


In terms of the velocity field and magnetic field we can determine 
the energy spectra; for the magnetic field, we find; 






( 19 ) 


and for the velocity field; 


<WS)I®> =5*;^<I|(S)|®> 







( 20 ) 


Similarly, the expectation value (B(k) • ;^-k)> is 


<B(k) • y(-k)> = f 


|a + 


X. . 


p® 

IJa 


-1 


= - |<iB(i')r"> 


( 21 ) 


Stunning over k gives the immediate consequence that (B • y) 

- .(p/a) <B^>. 

The conditions that the three expectation values of the 
invariants match their initial values. 
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e- 












-L 

k 


2 

k a 


7 - 


2 2 
g k 

iwt 


(22) 


are to be regarded as three simultaneous algebraic eqiiations which 
determine a, p, 7 in terms of e, P, A. Roughly speaking, 3/a measures 
the correlation energy of the magpietic field and velocity field. 

7/a is a more subtle parameter, but it will be seen presently that it 
meas'ires indirectly the ratio of long-wavelength magnetic energy to 
total energy to be expected for the equilibrium coefficients. 

A wide range of spectral behavior is possible for the different 
ranges of the initial parameters t, P, and A. One regime of impor- 
tance is the case P - 0. This is the case in which the cross correla- 
tion (B • y) = 0, either because B is initially uncorrelated with y 
or else is perpendicular to it. For this case, ji = 0 and the spectral 
densities simply reduce to 





( 21 *) 


P P 

Both (|B(k)j > ai«i <|;sKk)l ) approach constants ( equipartition 
spectra) at large jkj, but (!B(k)(^) can be peaked at |k| = 

o 

the longest wavelaigths, for 7 < 0 and a + 7 /k^^ smeJ.1. This 
regime always occurs if, for fixed ft and A, we allow to become 

large enough. This corresponds closely to the ’’Regime I” sharply- 
peaked energy spectra for the Navier-Stokes case (Kraichnan 1975 » 
Seyler et al* 1975)# in which large macroscopic vortices occur; here, 
it is ’’nagnetic islands” which appear. The vector potential spectral 
density, (|B(k)| /k >, is here dcaninated by the lowest values of 
I k) , while the extra energy is divided approximately equally among 
the highest k-values. The magnetic field configuration for this 
7 < 0 regime is dominated by the largest wavelength contributions, 
k = ( k ^^ ^ , 0 ) euid (0, k^^^ ) , and can be graphically described as 
’’magnetic islands" formed in the basic box, with essentially two 
regions of ciirrent densities of opposite sign dominating the flow 
pattern. From Eqs. (24) we can readily show that 


L <l|(k)|^> - E (|v()s)l^> 

k k 
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Since 2 k/ct is intrinsically positive^ it is a^^parent that the mag- 
netic field energy can have an expectati<m vidue which is greater or 
less than that of the kinetic energy, depending the eiga of 7. 

7 * 0 is a natural boundary in the parameter space of the problem. 
Fea* these (B • v) * 0 situations, the sl®a of 7 is basically deter- 
mined by the ratio A/ 1 , going from a sdnimiffl at A/g « to a 
maximum at k/t ® 0. 7-0 when 



and this ratio is, to a good approximation, k^^ ? for 

this value, (v > = (B ). For A/ 6 greater than this value (which -> 0 
as k^^ ■♦ *) , 7 < 0 and thus for fixed A and t, we always enter this 
regime if k „ is large enough. This regime is closely analogous 

niBA 

to Kraichnan's Regime I in the Navier-Stokes case. Note that Eq. 

( 24 ) predicts that when 7 < 0, (!B(k)|^> is greater than (lv(k)(‘ ) 
for all k, though only at the longer wavelengths is the difference 
significant. 

As 7 becomes more negative, the configuration increasingly 
resembles the limiting case of A/ g = which keeps all the energy 

locked into the two longest wavelength magnetic contributions. Two 
macroscopic electric current filaments of opposite senses occupy the 
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l 9 a«ic box In which the Foiffier analysis Is performed and have cross- 
sectional areas co^arable to it* To the extent that the periodic 
boundary conditions could be thouj^t of as simulating a finite physi- 
cal system, the external magnetic field to be observed would be that 
due to two parallel current carrying wlires with currents in c^sposite 
directions: a quadrupole field in two dimensicms. 

The limiting case of A/fc = 0 corresponds to 7 -* •*, a = [total 
number of allowed k's]/[total energy]. 'Ihere is no magnetic field 
energy at all, and a flat equlpartltion spectrum results for the 
velocity field. 

Between these two limits lies a continuous variation of con- 
figurations characterized, as 7 decreases, by a larger degree of 
Icsig-range order, with a more and more sharply peaked deposition of 
magnetic field energy at the longest wavelengths. The parallel with 
the formation of large vortices in the Navier-Stokes case, as the 
ratio of energy to entropl^ increases, is striking. 
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4. DISCUSSION 

Qualitatively, the phenomena predicted in the previous section 
ttre similar to a number of real situatlois in which a magnetic field 
of macroscopic dimensions has been observed to appear in a highly 
conducting fluid. Aouuig these, one may list the spontaneous genera- 
tion of megagauss ma^etic fields in laser -irradiated pellets, the 
persistence of "magnetic dynamo" fields produced in the earth’s core, 
and the ccnusolldation of current filaments In nun^rlcal simulations 
of energetic electron-beam situations (Lee and Lao^e 1913 )* For all 
of these situations, the geometry here assumed is an over-simplification, 
and detailed comparison with data would be premature, though the 
parallels with Lee and Lan^e, in particular, are interesting. 

Nevertheless, we find the emergence of the development of 
these spontaneously generated magnetic field structures from the 
mathanatics to have occurred in an unstrained and natural fashion. 

The phenanena would seem much less convincing to us if they had 
emerged as the end product of a time -dependent theory, which would 
undoubtedly have required any number of hard-to-Justlfy approximations, 
than they are, anerging as they have only frcxn arguments of an essen- 
tially thermodynamic character. At this stage, we are in the process 
of trying to discover the outlines of new dynamical effects, and the 
process of making them precise by including more realistic geon^tries 
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and finite dissipation terms Is likely to be as sloir as it has been 
for Havier-Stokes fluids. 

Hoirever^ we believe the i^enomena are interesting enough to 
stand on their own at this point as new predictions which nay be sub- 
ject to verificatim by numerical solution of the msffnetohydrodynamic 
eqwtims. What the t)»ory seens to be pointing toward is a class 
of macn^soopic iKsn-unifozm states of n<»-trlvial geometry, which the 
plaa&a regards as somehow natural and most probable. Because of the 
thermodynamic nature of their derivation, they almost certainly are 
stable. The class of c^iflgurations so described may weU. include some 
which various confinement schemes are and have been in search of by 
less systematic means. 
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